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ABSTRACT 

We analyze the effects of the detector response time on bolometric measurements of 
the anisotropy of the cosmic microwave background (CMB). We quantify the effect 
in terms of a single dimcnsionlcss parameter L defined as the ratio between the time 
the beam sweeps its own size and the bolometer response time. As L decreases below 
~ 2.5 the point source response of the experiment becomes elongated. We introduce 
a window function matrix based on the timestream data to assess the effects of the 
elongated beam. We find that the values of the window function matrix elements 
decrease slowly as a function L. Our analysis and results apply to other cases of 
beam asymmetry. For the High Frequency Instrument on board the Planck Surveyor 
satellite we show that for a broad range of L the ability of the experiment to extract 
the cosmological parameters is not degraded. Our analysis enhances the flexibility in 
tuning the design parameters of CMB anisotropy experiments. 

Key words: cosmic microwave background - methods: data analysis - methods: 
observational 



1 INTRODUCTION 

The anisotropy of the cosmic microwave background (CMB) 
radiation encodes a vast amount of information about struc- 
ture formation in the universe and about the values of 
the cosmological parameters (White, Scott & Silk f994). 
A number of groups are now making efforts to map the 
CMB at resolutions of f 0-20 arcmin. Two satellite missions, 
NASA's MAP and ESA's Planck Surveyor, are scheduled to 
be launched within the next decade. 

In a substantial fraction of ongoing and planned ex- 
periments the detector elements are bolometers. Bolometers 
are thermal detectors and thus have a response time that 
depe nds on a variet y of construction and operation parame- 
ters (Richards 1994). Time constants of bolometers currently 
employed on CMB experiments are on the order of 10 msec. 
The baseline design for Planck's High Frequency Instrument 
(HFI) bolometers is between 2 and 5 msec. 

The bolometer response time puts a practical constraint 
on the speed with which a bolometric experiment can scan 
the sky. Scanning too fast decreases the sensitivity of the 
experiment to structures near the beam resolution. On the 
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other hand, the scan speed and the detector's 1/f noise knee 
determine the number of time-ordered pixels over which the 
noise is uncorrelated. A faster scan speed allows a more 
robust characterization of the instrumental noise and bet- 
t er recovery of t he sky signal o ver a larger angular range 
( [Tegmark 19974 |Tegmark 1997b| ) . Thus there are important 
trade-offs between the detector noise level, time constant, 
beam size and scan speed that must be considered when 
optimizing the experiment. In the case of the Planck Sur- 
veyor instruments, for example, a balance has to be found 
between the requirements of the low frequency HEMT-based 
instrument (LFI) and the bolometric high frequency instru- 
ment. The LFI prefers a fast scan speed because both its 
detector response speed and 1/f noise knee are higher than 
for the HFI. However, if the scan speed is too fast HFI's 
performance at high resolution will degrade. 

The window function matrix, and its trace, the window 
function, quantify the sensitivity of a CMB anisotropy ex- 
periment as a function of angular scale. The full matrix con- 
tains information about all possible pixel pairs; the simpler 
window function, We, reduces this to the RMS anisotropy 
probed by the experiment. A number of authors have ana- 
lyzed how the beam size and scan st rateg y determine the 
window function (see White et al. (1994) and references 
therein). It has traditionally been assumed that the bolome- 
ter response time can be neglected. With a window function 
at hand, algorithms have been developed to asses the accu- 
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racy with which a given experiment can extract the cosmo- 
logical parameters (Jungman et al. 1995; Bond, Efstathiou 
& Tegmark 1997; Zaldarriaga, Spergel & Seljak 1997). All 
of these calculations assume that the experiment illuminates 
the sky with a symmetric beam and that the noise in the 
experiment is uncorrelated, and hence that off-diagonal ele- 
ments of the window function matrix are unimportant. 

In this paper we analyze how the window function ma- 
trix of a bolometric experiment depends on the combina- 
tion of bolometer response time, beam size, and scan speed. 
When the detector response time cannot be neglected, the 
effective shape of the beam becomes elongated and the form 
of the window function matrix becomes a complicated func- 
tion of the scan strategy and pixelization. Furthermore, the 
standard techniques used to assess the accuracy with which a 
given experiment can extract cosmological parameters from 
the data are no longer strictly valid since they assume that 
the beam is symmetric. We overcome the difficulties asso- 
ciated with constructing a window function matrix and as- 
sessing the performance of the experiment in the presence of 
an asymmetric beam by introducing a window function con- 
structed from data points pixelized in the time domain. We 
use this 'temporal' window function matrix to evaluate the 
window function as a function of the detector response time 
and to estimate the magnitude of the necessary corrections 
to the off-diagonal elements. The discussion focuses on the 
potential reduction in the window function matrix response 
near the nominal resolution of the experiment. Our analy- 
sis and results will be applicable to any asymmetric beam 
elongation that decreases the nominal resolution in one di- 
rection. As a special application, we calculate the ability of 
Planck-HFI to extract cosmological information as a func- 
tion of different values of the instrumental parameters. 

In Section^ we derive the effect of the detector response 
time on the point source response of the experiment. We 
then derive the necess ary changes to the zero lag window 
function (Sectio ns |3.l| , 3.2) and asse ss Planck-HFI's perfor- 
mance (Section |3.3[ ). In Section 3.4 we concentrate on the 
full window function matrix in the presence of the beam 
asymmetry introduced by the detector response time. We 
discuss and summarize the results in Section W. 



2 BOLOMETER TIME CONSTANT AND 
POINT SOURCE RESPONSE 

A bolometer is a thermal detector in which absorbed radia- 
tion is converted to heat, causing a temperature change that 
is proportional to the absorbed energy. A thermistor is used 
to measure the temperature change. The temporal response 
of bolometers is an exponential; if the bolometer tempera- 
ture is To at time t < 0, then upon a step function increase 
in input power at time t = 0, 



(2) 



T(t > 0) = T + AT(1 - exp(-t/r)). 



(1) 



The time constant r depends on the heat capacity of the 
absorbing medium, on the detector's thermal conductivity to 
i ts mounting st ructure, and on properties of the thermistor 
( Richards 1994 ) . The exponential response time r is a single- 
pole low-pass filter in the frequency response of a bolometric 
experiment. The filter is given by 



1 + (ujt) 2 ' 

The amplitude response of the filter has a —3 dB pointj^j at 
a frequency /bolo = l/(2nr). 

In an observing strategy where the optical beam is 
scanned across the sky the bolometer time constant can be 
viewed as a single-pole low pass filter in the spatial domain. 
Since all of the effects we will be discussing are relevant only 
at small angular scales, we approximate the region of interest 
on the sky as flat and use Fourier transforms on the plane. 
If we write u> = k ■ v, where k is the spatial wavenumber 
being probed, and v is the angular velocity of the pointing 
on the sky, we get for the response in k-space 

1 — i kvr cos 6kv 



F(k) = 



1 + (kvr COsdkv) 2 



(3) 



where 6kv is the angle between the vectors k and v. F is a 
complex filter. Its amplitude response is 



\F(k)\ = 



1 



+ (kvr) 2 cos 2 9 k 



and its phase is 

(f> — tan _1 ( — kvr cos 



(4) 



(5) 



The amplitude response of the filter has a —3 dB point at 
k^dB = 1/(vt cos 9kv) which depends on the instantaneous 
scan speed v. From now on we will implicitly assume that, 
where relevant, the scan speed is constant. The approxima- 
tion of a flat sky is valid for regions smaller than ~ 20 x 20 
deg 2 . Within this approximation k — £, where £ is the mul- 
tipole number of the Yi m spherical harmonics. 

Assume that the optical system of a bolometric experi- 
ment produces a symmetric beam B(r, a), where a is a mea- 
sure of the size of the beam (e.g., the width of a Gaussian 
beam), and r locates it. The position r is defined relative 
to coordinates located on the region of interest on the sky. 
If the beam is swept across the sky at a constant velocity 
v the point source response is the convolution of the beam 
and the temporal response of the detector: 



5(r,v,r, a) = 
equivalently, 



B( r — r' , <x) F(r, v, r) d 2 r'; 



S(r, v, r,a)= / exp(-ik • r) B(k, cr) F(k, v, t) 



ri 2 k 



(6) 



(7) 



where, to simplify notation, we use B(k) and F(k) as the 
Fourier transforms of B(r) and F(r). Since F(k) is complex, 
the resulting response S is both attenuated and phase shifted 
compared to the case r = 0. A dimensionless figure of merit 
that quantifies the level of attenuation and phase shift is L, 
defined as 



VT r 



O/V _ /bolo 



fm 



(8) 



Since a/v is the time it takes the beam to cross a distance 
equal to its width, L>1 corresponds to the case where the 
time constant is short compared to this crossing time and 
spatial filtering is negligible. L < 1 corresponds to the case 

§ dB = 201og 10 (|F|) 
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Figure 1. Top panel: The point source response of a Gaussian 
beam as a function of L = a/(vr), where <r is the width the 
beam, v is the scan speed, and r is the bolometer time constant. 
For t > 1 there is little change in the shape of the beam. For 
L < 2 the point source response is phase shifted and attenuated. 
Bottom panel: For values of L > 2 the effective shape of the beam 
in the direction of the scan is well approximated by a Gaussian 
with a width a 2 d ~ a 2 + (vt) 2 . Shown are the effective shape 
(full curve) and the Gaussian fit (dash-dot) for L = 2. 



where the time constant is much slower than the crossing 
time and spatial filtering is expected to be significant. In 
the frequency (1/time) domain, L is the ratio between the 
bolometer time constant low-pass filter /boio and the width 
of the point s ource response with r = 0, / snm = v/(2ira) 
(Hristov 1994). For a Gaussian beam / SC an is at the 1/e 
point of the point source response. 

The upper panel in Figure [j] shows the point source re- 
sponse of an initially Gaussian beam for a range of L values. 
For a nominal two-dimensional Gaussian beam, 



B(r,a) = 



2tt(J 2 



-r 2 /(2cx 2 ) 



(9) 



(10) 



As expected, small values of L result in a larger effective 
beam in the direction of the scan. The shape in a direction 
perpendicular to the scan direction is not affected. Overall, 
the larger effective beam size will result in a degraded sen- 
sitivity to spatial structures which are on the scale of the 
nominal resolution. 

For values of L Jj 2 the spatial filter (Eq. ^) can be ap- 
proximated as a Gaussian. When combined with a nominal 
Gaussian beam of width a the Fourier transform, Eq. [j], can 
be solved exactly, and one finds that the effective beam in 
the direction of the scan has acquired an additional Gaus- 
sian of width a 2 olo ~ (vt) 2 . The total beam width in the 
direction of the scan is now 



2 

&ld 



a 2 + (vt) 2 = cr 2 (l + L~ 



0.4 



(11) 

The lower panel of Fig. |l| shows the effective beam shape 
and the Gaussian that approximates it for L — 2. The re- 
sulting 2-dimensional beam is elliptical with a width aid in 
the direction of the scan and a in a direction perpendicular 
to the scan. 



3 WINDOW FUNCTION AND 

COSMOLOGICAL PARAMETER 
ESTIMATION 

We now quantify the effect of the enlarged asymmetric beam 
on the window function of the experiment and its effect on 
the estimation of cosmological parameters. We ignore the 
spatial phase shifts of the signal. In most real situations the 
phase shift will either be small compared to the beam size 
or the need to accurately reconstruct the beam location (for 
e.g., making an accurate map) will force experimenters to 
determine and correct for it. 



3.1 Temporal Window Function 

We begin with the temporal data stream which consists of 
data labeled by both position on the sky and time. The 
signal at a given time t = 1 . . . Nt is a convolution of the 
beam pattern with the underlying sky signal: 



Si 



AT 

d 2 x -F t (x) — (x) 



aim 



(12) 



where we have transformed to spherical harmonics in the 
second equality. The kernel T encodes the beam shape and 
its location on the sky at time t. If we consider small areas 
of sky, on which the time-constant effects will be relevant, 
we can again consider the sky to be flat and approximate 
the spherical harmonics with Fourier transforms. Then the 
signal is 



Si 



d 2 k 



-a(k)^(k), 



(2tt) 2 
with 

(a(k)a(k')) = (27r) 2 <5 2 (k + k')C(fc), 
and 



(13) 



(14) 
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C e ~C(k)\ e ^ k . (15) 
The temporal signal correlation matrix is 

k 2 C{k) 



S u > = (atsf) = / dlnkW tt '(k) 



2tt 



(16) 



with the temporal window function matrix 



Wu>{k) = J -^e- tKr "' COSHk U tt ,{^)B\k) (17) 

where we have split off the effect of a symmetric beam B(k) 
and other asymmetric parts U. Here, r tt < is the angular dis- 
tance between pixels observed at times t, t' . An experiment 
with U — 1, (i.e., a symmetric beam, and no chopping), has 

W tt '(k) = Mkr tt ,)B 2 (k). (18) 

For t = t' we get the zero-lag window function 

Wu(k) = B 2 (k). (19) 

With the spatial asymmetry of the beam, induced by the 
bolometer time constant (Eq. Q), the temporal window func- 
tion matrix becomes 

W tt ,{k) = B 2 (k) ( d^e-*"-*™*"* 



x [i + (kvr) 2 cos 2 (e k -e t )] 1/2 

x [l + (fc^r) 2 cos 2 (6» fe -0 t /)]" 1/2 , 



(20) 



where 9± , 9 t i are the angles of the velocity vector v at times 
t and t' . 



3.2 Zero Lag Window Function 

To obtain the traditional window function of the experiment, 
which encodes contributions only to the RMS signal, we set 
t = t' in Eq. pel, giving 



Wu(k) 



B 2 ( k ) I t 1 + (kvr) 2 cos 2 (0 k -0 t )] 1 



B 2 (k)- 



1 



where 



y/l + {k/k r ) 2 

^ = J_ _ 1.8 ■ 10 5 /lms\ ( ldeg/s 

VT 7T V T / 



V V 



(21) 



(22) 



The original window function, which was composed of the 
beam spatial filter, W T =o(k) = B 2 (k), is now multiplied by 
an additional single-pole low-pass filter, 



R(£) ~ R(k) = W(k)/W T=0 (k) 

= [1 + (k/k T ) 2 ]~ 1/2 ~[1+ (£/l T ) 2 ]- 1/2 . 



(23) 



The —3 dB point of R(£) is at £ = l T . The beam filter has a 
low-pass cut-off at l a = 1/a. The shape of the filter depends 
on the shape of the beam. ( A Gaussian beam with width 
a has an exponential low pass B 2 (k) = exp (— £ 2 a 2 ).) Thus 
the effects of beam elongation will be important when 



e T < £ a => l = -f- <, i 



(24) 



Figure |2j shows the zero lag window function for an originally 
Gaussian symmetric beam and for various values of L. 



.00 



0.10 



0.01 













- 








- 


- L (fro 
: L > 1 

L = TT 

L = 2 
L = 1 
L = 


n right to 
30 

.5 


eft) \ \ \ 







0.0 



0.5 



1.0 1.5 
I (units of O 



2.0 



2.5 



Figure 2. The zero lag window function, Eq. ^l], with an initially 
Gaussian beam for a range of L values. The horizontal axis is in 
units of £(j = X/cr. 



We point out that the window function calculated from 
the timestream, Eq. ^l], is the same as the one calculated 
from a map, 

Wt = j r y £ i W„{£), (25) 

P N p 

where p denotes a pixel location on the sky. The window 
function, in any representation, just measures the sensitivity 
to the RMS temperature, and does not take into account 
pixel-pixel correlations. 

For a reasonable range of v and r, £ T ^ 800, below which 
most of the cosmological information resides in CDM-like 
models. For these cases, and for £ £ T , the single-pole filter 
Re can be usefully approximated by a Gaussian. By Taylor 
expanding (1 + £ 2 / 1 £ 2 )~ 1 ^ 2 and matching with the first term 
of a similar expansion for a Gaussian we find that 

-e 2 o-' 2 _ -e 2 (<7 2 +o-' 2 ) 



W g (£)~B 2 (£)e- 

where 

a = J_ _ (w) 2 
a ~ 2£ 2 2 ' 



(26) 



(27) 



and we have specialized to a Gaussian beam in the equality. 
This can be thought of as approximating the nearly ellip- 
soidal effective beam by a circular Gaussian beam with a 
slightly larger width. This new width, 

(28) 



2 

"2d 



2 (vt) 
(J' + - 



:a 2 (l + L- 2 /2) 



is larger than the original width a and smaller than the 
one dimensional effective width that we found in the pre- 
vious section; see Eq. [H]. That is, we find a degradation in 
the window function response compared to the original a 2 
width, but not as severe as indicated by the beam elonga- 
tion in the scan direction, a\ d = a- 2 (l+L -2 ). This is because 
beam elongation occurs primarily in one direction whereas 
the window function encodes information from all spatial 
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directions. Generally, as can also be observed by Comparing 
Figs. and | the effects of small values of L are stronger in 
the spatial domain than they are on the window function. 

The effective Gaussian filter exp(— £ 2 (a 2 + a' 2 )) is the 
geometric mean of two Gaussian filters corresponding to two 
symmetric beams (White 1997), 



B 2 (k) 



^/B 2 (k)B 2 ld (k). 



(29) 



One has a width aid corresponding to the effective beam in 
the scan direction and the other has a width a corresponding 
to the nominal beam size. The Gaussian approximation is 
good to within 10% for all £ less than .88 l T . As an example, 
l T = 1680 for t — 5 msec, and v — 6 deg./sec. which are 
the time constant for Planck's 100 GHz bolometric channel, 
and the satellite rotation rate, respectively, and thus the 
approximation is valid for £ < 1480. 



3.3 Estimation of Cosmological Parameters 

Given the window function W(k) (Eq. ^l|), we use the appa- 
ratus that^jM£l2^er^ (1996) and 
others (Bond et al. 1997; Zaldarriaga et al. 1997) to estimate 



its effect on the determination of cosmological parameters. 
We specialize the discussion to the experimental configura- 
tion of the High Frequency Instrument (HFI) on board the 
Planck Surveyor satellite. 

Each of the panels of Table [j] represents a different cos- 
mological model. For each of the models we calculate the 
error in the determination of a variety of cosmological pa- 
rameters for three values of the parameter L. Lo = 2.67 is 
chosen to represent the value presently designed for the HFI 
instrument. 

Within the bounds of validity of the analysis (see be- 
low), we find that for the case of the HFI choosing L val- 
ues which are a factor of 2 or 5 smaller than the nomi- 
nal Lo would not significantly affect the ability to extract 
the values of the cosmological parameters. Note that the 
effect of the smearing is almost completely negligible for 
L ?Z Lo/2 ~ 1.33 and begins to be noticeable with the three 
significant figures in the table at L ~ Lo/5 ~ 0.5. We defer 
a more detailed discussion of this result to Section ^. 

The results in Table [j] should be read with some cau- 
tion. The forecasts of parameter errors assume a form for 
the cova rianc e matrix (more precisely, the Fisher Matrix; see 
Section fc.4.l| ) which is appropriate for an experiment with 
a symmetric beam and uniform, uncorrelated pixel errors. 
Since in our case the beam is manifestly not symmetric we 
expect these forecasts to be only approximate. Furthermore, 
our analysis considers quantities in the time domain. The re- 
lationship between pixel- and time-domain is not completely 
straightforward in the presence of beam asymmetry; this is 
a topic of current research. In the limit of many crossings 
through the same pixel in many different directions, the re- 
quirement for an effectively symmetric beam is satisfied, if 
we consider quantities in the pixel, rather than time, do- 
main. In other cases it will not be satisfied and care mu st 
be taken in interpreting the results in Table |l|. White ( 199 



has noted that it may be possible to get a pixel-domain win- 
dow function that is narrower than the time stream window 
function that we use here by certain signal-to-noise weight- 
ing techniques. More generally, it has been noted that the 
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Table 1. The three columns of numbers quantify the accuracy 
with which the Planck Surveyor HFI will be able to determine 
a variety of cosmological parameters as a function of L = a/vr, 
where a is the width of a Gaussian beam, v is the scan speed, 
and r is the detector time constant. The parameter Lo = 2.67 
is appropriate for the baseline design of the HFI. The calcula- 
tion assumes 4 frequency channels at 100, 143, 217, 353 GHz, 
Gaussian beams with FWHM 10.6, 7.4, 4.9, 4.5 arcmin, nominal 
sensitivities of 4.9, 4.8, 9.7, 40/iK /pixel., respectively, and a scan 
rate of 6 deg./sec. The notation 'rel' is for fractional errors with 
respect to the model. Due to a degeneracy between f! V ac and Qq, 
we hold the curvature, 1 — f!o — Slvac — — f2„, fixed for the 
error determination, except in the determination of Qa itself, for 
which we fix f2 V ac- 



least-squares method of Wright (1996), which is a variant 
of the method originally used to make the C OBE/DMR 
maps (see also (Tegmark 1997a; Tegmark 1997b)), produces 
a map containing all of the sky information of the experi- 
ment. The method, as implemented in these references, as- 
sumes a known timestream noise matrix N tt / and uses it to 
estimate the temperature at a pixel p given measurements 
of the temperature at p taken at different times. With an 
asymmetric beam, individual measurements taken at p no 
longer measure quite the same quantity due to the differing 
beam orientations. This can be accounted for at the cost 
of greatly complicating the matrices used in the solution to 
the problem. Such a treatment is beyond the scope of this 
paper. For the results of Table [j] we used timestream quan- 
tities and the RMS window function W(k). We will return 
to the question of beam asymmetry and assess its impact on 
the window function matrix, and on the ability to extract 
cosmological parameters, in the following section. 

Finally, the standard analysis, which leads to the results 
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in Table [l], makes several additional assumptions that need 
to be highlighted. The actual likelihood as a function of the 
cosmological parameters will be highly non-Gaussian, but 
the calculations assume Gaussianity (at least near the peak 
of the distribution). In addition, the results assume that 
the removal of foregrounds and systematic effects proceeds 
without appreciable degradation of the satellite's sensitivity. 



write the Fisher matrix in terms of symmetric-beam quan- 
tities: 

1 



Ft 



2 Tr 

rer e > 
o? 



(aS + N)- L r t W {aS + N)~ x r t >W 



(s + N/ a y 1 w\s + N/ay^w 1 ' 



■F u ,(S,N/a) 



(35) 



3.4 Full Window Function Matrix 

The analysis in terms of the zero-lag window function is 
an approximation that neglects pixel-pixel, or beam-beam 
correlations. We now assess whether sweeping an effectively 
asymmetric beam across the sky significantly affects this ap- 
proximation. We first parameterize the relative power spec- 
trum errors 8Ct/Ce in terms of the symmetric beam errors 
and a. The parameter a quantifies the 'average' ratio in the 
beam-beam correlations between the symmetric and asym- 
metric beam cases. We then assess the magnitude of a. 



3.4-1 Off-Diagonal Elements and Parameter Estimation 

The (ensemble average) error on the cosmological parame- 
ters bi is given by 



(5k 5b i ,)=FZ,\ 

where F it i is the Fisher Matrix given by 



iV- 2 TrlG C 



(30) 



(31) 



Here C = C pp > is the total correlation matrix (in pixel space, 
although the arguments hold for the timestream data as 
well) of the data in question, which has contributions from 
both signal and noise: C = C(S, N) = S + N. In particular, 
the signal part of the matrix is given by 



(32) 



Thus, if we parameterize by the power spectrum itself the 
derivatives are just 

dC _ OS _ 21 + 1 , 

dCt ~ 



dC e 



-Wt , 

4tt pp 



(33) 



Here we write the full window function matrix as W, 



pp' 



and the usual 'window function' as We = (l/iVp)Tr(W ) = 
(1/Np) Wp p . The Fisher matrix for any other set of pa- 
rameters can be recovered from that of the power spectrum 
by the appropriate Jacobian matrix: 



E 



dCt 
db % 



' dbi, 



(34) 



That is, we can write the Fisher matrix as approximately a 
factor times the Fisher matrix for a symmetric beam, with 
the noise variance degraded by 1/a. We might expect this 
factor, r^ry/a 2 , to be of order one. 

In particular, for an experiment with uniform noise, we 
can write the Fisher matrix for an experiment observing 
a fraction / s k y of the sky with a symmetric beam Be as 
(approximately) 

-I -2 



f u < ^ (e + i/2)f sky 



C e + 



wB 2 



(36) 



Here, w = Tr(AT _1 ) is the total weight of the experiment, so 
the effect of the beam asymmetry is to make the replacement 
w — > aw; when the noise term dominates, the error on Ce is 
increased by 1/a. 

3-4-2 Off-Diagonal Elements of the Window Function 
Matrix 

We now assess the magnitude of the off-diagonal elements 
of the window function matrix in the presence of beam elon- 
gation due to small values of the parameter L. We will be 
interested in the ratio of Eq. |2(j, which is the general expres- 
sion for the window function with an asymmetric beam, to 
Eq. which assumes a symmetric beam 

Rtt'(k) = W tt ,/{B 2 (k)J (kr w )} 

1 



J (kr t 



ddk -ikr t 
27T 



x [l + (kvT) 2 cos 2 (e k -0t)]~ 
x [l + (br) 2 cos 2 (0 k - Of)] 



'1/2 
~l/2 



(37) 

Note that the original symmetric beam B(k) drops out of 
this expression. When this ratio is small the beam elongation 
is having a large effect. While this ratio cannot be calculated 
analytically, we can gain much insight into the time-time 
correlations by considering the range in which kvr = i/l T 
is less than one. We then compute the integral numerically 
for larger values of kvr. 

Expanding Eq. ^ in terms of a = (kvr) 2 , redefining 
the zero point of the angular integral, and working to order 
a 2 we find 



Now we can determine the effect of a window function 
with an asymmetric beam on parameter errors. We refer all 
the requisite quantities to those with a symmetric beam. 
Define WjLi ~ reW pp r, where W pp > is the window function 
for a symmetric beam and re ~ R(£) is an appropriate pixel 
average of the ratio R we defined in Eq. |2^. Also define 
S pp i — aS pp i where again S is the quantity for a symmetric 
beam, and a is thus a weighted average of the re. Since the 
form of W will not be the same as that of W, these ratios 
are of course only approximations. With these definitions we 



Rtt>(k) 



+ 



Jo(kr tt > 
% [ Sin2( 



d0 ikr s 

— e 

2tt 



6)+sm 2 {6 t , -0)] 



! [lsm\0 t -0) 

2 1 • 2 rn 
+ a - sin f t 
4 

+ |sin 4 (^- 



lsin 2 (6» t / 



•«o] 



(38) 
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By contour integration, this can be broken up into a sum of 
spherical Bessel functions: 



Ru'Joikr) = Jo(fcr) 




cos 4 Ot + cos 2 9t cos 2 8 t i + ^ cos 4 ( 



x [3J (fcr) - 4J 2 (fcr) + J 4 (fcr)] 



(9 cos 2 Ot sin 2 # t + cos 2 t sin 2 (9 t / 
+4 cos 8t sin St cos Q t ' sin #f 
+ sin 2 St cos 2 #f + 9 cos 2 t i sin 2 # t ' 
x [Jo(fcr)- J 4 (fcr)] 

3 . 4 a . • 2 fl - 2/j .3 .4 

— sm t?t + sin Ot sm t? t / + — sin ( 





Case 1 






Case 2 



Case 3 



Figure 3. Three beam configurations that contribute to the off- 
diagonal elements of the window function matrix. The configura- 
tions shown span the range of interesting topologies. 



x [3 Jo (Ax) + 4J 2 (fcr) + J 4 (Ax)] 



(39) 



We note that for tr>l 



Jl(x) > - ! 

X 

and thus Eq. 
Rtt< = 1 



(-t) 



~2~ 

becomes 



(40) 



- [2 + cos(26> t 



+ cos(26l t /)] 



+ — (3 cos 4 I 



+ 2 cos 2 6t cos 2 ( 



+ 3 cos 4 6> t ,)(.41) 



We examine Eqs. [39] and ^ for three different configu- 
rations of asymmetric beams that are the interesting cases 
in the range of possible topologies. Fig. ^ shows the three 
configurations. In Case 1 the beams are separated along the 
direction of the scan. In Case 2 the beams are separated 
perpendicular to the scan, and in Case 3 the separation be- 
tween the beams is parallel to the scan direction for one 
beam, and perpendicular for the other beam. For each one 
of these configurations we fix the distance scale r at two val- 
ues, and plot in Figure ^ the behavior of R tt ' as a function 
of fc. We let k vary from k = up to fc ma x = l/(2«r) so that 
the approximation a < 1 is valid throughout. 

Case 1, shown in the upper left panel of Figure ^, is ob- 
tained by taking Ot — 9 t i = in Eq. KEN. The periodic spikes 
in the value of R t t' are due to the slight shift between the 
zeros of Eq. p9| and Jo(fcr). Ignoring these spikes we observe 
the following properties of R t t> ■ Near k = 0, i.e., for very 
long wavelengths, there is no difference whether the beams 
are elongated or not. As k increases R t t' decreases for both 
values of r indicating a smaller off-diagonal elements com- 
pared to the symmetric beam case. There is less attenuation 
for small r as compared to large r. 

To understand this we note that only those k vectors 
that are directed near the line between the two beams con- 
tribute to the angular average performed in Eq. Since 
the asymmetric beams are most widened in that direction 
large fc's are effectively smeared (hence R t t' < 1)- Larger 
attenuation is obtained at large beam separation because 
a narrower range of wave vectors, and only those that are 
aligned in the beam elongation direction, contribute to the 



1.1 
1 

, 09 
0.8 

0. 



Case 1 



m 



I k 
I I 

I I- 



1.1 

1 
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X 

0.8 
0.7 
0.6 
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Case 2 
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0.2 0.4 0.6 0.8 1 
k / k m „ 



Case 1 (Numerical) 




0.2 0.4 0.6 0.8 1 

k / k. 



Figure 4. R u >(k) = W t v I 'B^(k)Jo(kr) for the three different 
topologies illustrated in Fig. bl The upper left panel is Case 1, 
the upper right panel is Case 2, and the lower left panel is Case 
3. In each panel r = l/fcmax (dotted) and r = 20/fc max (long 
dashed), where fc max = l/(2vr). The periodic spikes are due to 
the shift between the zeros of Eq. ^9] and Jo(fcr), and they can 
be ignored. Note that in these plates the vertical axis starts at 
R tt / = 0.6. In the lower right panel we show Case 1 again for a 
larger span in k and for r = 20/fejpax. These results were obtained 
by numerically integrating Eq. B7J. In the range of overlap in k the 
top left and bottom right panels agree. 



angular average. Within this range of fc's the largest attenu- 
ation compared to the symmetric beam case is ^ 20%. This 
limiting value is also obtained from Eq. |H] with 8 t = Of = 
which gives R t t> ~ 1 ~ a + a 2 — 0.8. 

Case 2 is obtained by taking 0t = t > = 7r/2 in Eq. |39[ 
The upper right panel of Fig. ^ shows R tt '(k) for this case. 
Similar to Case 1, when fc = we obtain the symmetric beam 
case, -Rtt' = 1. For r = 20/fc m ax there is no change in R as 
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k increases, and there is ~ 10% change for the r = l/fc max 
case. 

The physical explanation is similar to the reasoning in 
Case 1. The beams are narrowest in the direction of their 
separation. For large r only those k vectors that are aligned 
with the narrow dimension of the beam contribute to the an- 
gular average in Eq. Hence there is effectively no differ- 
ence between the symmetric beam case and the asymmetric 
case. The large r limit as obtained from Eq. ^] is R tt i = 1. 
At small r values k vectors which are in the direction of the 
beam elongation contribute to the angular average. In this 
case R tt i is smaller compared to either the symmetric case 
or to the large r value. 

Case 3 is obtained by taking 0t — and 9 t i — n/2 in 
Eq. The lower left panel of Fig. ^ shows R u i (k) for this 
case. For arguments similar to the ones in the previous cases 
we expect, and indeed observe, little difference between the 
two values of r. We also observe an overall attenuation at 
fcmax of ~ 10%. This agrees with the asymptotic value for 
large r, R tt ,{k) « 1 - a/2 + 3/8a 2 = 0.9. 

Within the approximation kvr — £/£ T < 0.5 the largest 
reduction in the off-diagonal elements occurred for case 1, 
and for large values of r. For this case we continue the anal- 
ysis for larger values of kvr by numerically integrating Eq. 
|37| . The result of this numerical integration is shown in the 
lower right panel of Fig. ^. We find that the analytic ap- 
proximation is excellent for kvr = £/£ T < 0.5 and that when 
£/l T ~ 1 the magnitude of these particular off-diagonal ele- 
ments is about 50% of their values in the symmetric beam 
case. 

We emphasize that we are considering here the effect of 
the asymmetric beam relative to the symmetric beam case. 
In Case 1 the asymmetric beam causes the greatest relative 
effect at large separation. However, the absolute contribution 
of these matrix elements to the signal at large £ will be small 
because large separations sample mainly small £ values (that 
is, the matrix elements are suppressed for fcr ~ fr > 1). 
We thus find that for the £ range where beam elongation is 
important, i.e., large £, the matrix elements are not atten- 
uated significantly compared to the symmetric beam case, 
certainly less than the 50% reduction observed for the most 
seve re cas e at £/£ T — 1. Based on the arguments of Sec- 
tion 3.4.1 we expect a similar increase in the relative errors 



SCe/Ce- The change in the relative error expected for any 
one of the cosmological parameters will be similar or smaller, 
and thus the major conclusions of Table |l| will not change. 



factor of two smaller than that derived from the point source 
response (Eq. 

Our analysis of the performance of Planck-HFI indicates 
that choosing L ~ 0.5 will not be detrimental to the ability 
to extract cosmological parameters. This is a consequence of 
HFI's small beam size and high sensitivity. The HFI, with 
its nominal £ a = 1/a which range from 725 to 1800, will 
determine the CMB power spectrum with high accuracy over 
a very broad range of £. A degradation of its performance at 
high £ values does not diminish its capability significantly. It 
will generally be the case that the ability of any experiment 
with £cr 800 to extract cosmological parameters will not 
be affected significantly with L J; 0.5, if the universe is any 
variant of the CDM model. 

There are strong arguments for maintaining high resolu- 
tion capability for a CMB experiment. We do not know that 
the universe is any variant of the cold dark matter model, or 
indeed that the fluctuations are Gaussian. High resolution 
will undoubtedly be important to test this latter assump- 
tion. Also, the ability to remove foregrounds and to detect 
point sources relies to some extent on high £ information and 
the ability to observe small-scale non-Gaussian structures. 

If we are interested in achieving the maximum possible 
resolution, we wish to minimize a quantity like the effective 
2-d beamwidth a\ d ~ cr 2 (l + L~ 2 /2). On the other hand, 
other considerations will make a smaller L more desirable. 
A fast scan speed increases the scale over which instrument 
noise is uncorrelated. Such broad temporal and/or angular 
scales can be used to better characterize the instrument's 
noise, to increase the £ range to which the experiment is 
sensitive, or for a bet ter recovery of t he underlying two di- 
mensional sky signal (Tegmark 1997a). The form of 02d ar- 
gues that we quickly reach a regime of diminishing returns 
as L > 1. For a given experiment the final determination 
of L depends on a global optimization which includes the 
expected instrument noise characteristics and the targeted 
science. Our analysis indicates that 1 ^ L ^ 2.5 provides 
a range over which the small angular scale performance is 
not affected significantly, while allowing flexibility in tuning 
either the bolometer response time and/or the scan speed. 
For the HFI in particular, relaxing the value by a factor of 
~ 2, to L ~ 1.5, may increase the dynamic range over which 
the instrument and satellite parameters can be varied. Some 
experiments may choose to use even smaller values of L. 



4 DISCUSSION 

Our point source response analysis shows that using a com- 
bination of bolometer time constant, scan speed, and beam 
size such that L — a/vr = £ T /£cr J5, 2.5 results in relatively 
small increase in the effective beam size, less than 16%, and 
the spatial phase shift is less than 0.2 a. Because the win- 
dow function encodes information from the nominal beam 
dimension as well as from the enlarged one the effect of a 
smaller L is less pronounced on the window function than 
it is on the point source response in the scan direction. In 
the Gaussian approximation the fractional effective increase 
in beam size derived from the window function analysis is a 
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